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In [1] we introduced a class of positive recurrent Markov chains, named 
tame chains. A perfect simulation algorithm, based on the method of dom- 
inated CFTP, was then shown to exist in principle for such chains. The 
construction of a suitable dominating process was flawed, in that it relied 
on an incorrectly stated lemma ([1], Lemma 6). This claimed that a geomet- 
rically ergodic chain, subsampled at a stopping time a, satisfies a geometric 
Foster-Lyapunov drift condition with coefficients not depending on a. This 
is true if a is a stopping time independent of the chain, but not if this 
independence does not hold. Reference [1], Lemma 6 is therefore false as 
stated. 

We now indicate a corrected construction of a dominating process. As 
described in [1], Section 3.1, the process D is defined by starting with a 
process Y and pausing it using a function S. In the following modified con- 
struction this is simplified by taking S = F, where F is the function taming 
X. We restate [1], Theorem 16, and give a shorter proof, which avoids the 
faulty Lemma 6 but pays a price in terms of consequences for the perfect 
simulation algorithm of Section 3.3. The discussion of tameness (Section 4) 
is unaffected. 

Theorem 16. Suppose X satisfies the weak drift condition PV < V + 
61c, and that X is tamed with respect to V by the function 

\Xz 5 ], z>d', 
1, z<d', 



F{z) 



with the resulting subsampled chain X' satisfying a drift condition PV < 
(3V + t/lty<d']> with l°g/3 < <5 _1 k>g(l — 5). Then there exists a stationary 
ergodic process D which dominates V(X) at the times {o~ n } when D moves. 
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Proof. Suppose that D an = z, and that V(X CTn ) = V(x) < z. We wish 
to show that D an+1 can dominate V(X an+1 ), where cr n +i =ff n + F(z) is the 
time at which D next moves. Domination at successive times <7j at which D 
moves then follows inductively. For simplicity in the calculations below we 
set a n = 0. 

First choose (3* > f3 such that 

(1) log ,3 < log /T <^ 1 log(l-5). 

Our aim is to control Ex[^(^f(,z))]> recalling that F(z) is deterministic and 
that F(V(x))<F(z): 

E x [V(X F{z) )}=E x [V(X F{v{x)) )}+E x [V(X F{z) ) - V(X F{V{X)) )} 
= E X [V(X[)} +E x [V(X Hz) ) - V{X F{v{x)) )} 
< (3V{x) + b'l [v(x) < dl] + b[F(z) - F(V(x))} 

<f3z + b' + b(X + l)z s 

(2) <(3*z forz>h*, 

where h* < oo is a constant chosen sufficiently large for inequality (2) to 
hold. The first inequality in this sequence holds due to the drift conditions 
satisfied by X' and X. The second follows from the definition of F and the 
assumption that V(x) < z. 

Now define the process Y = h* exp(U), where U is the system workload of 
a D/M/l queue with arrivals every log(l//3*) time units and service times 
being independent and of unit Exponential distribution. As in the original 
proof of Theorem 16, Y may be paused using F to obtain the process D 
which is positive recurrent and has a proper equilibrium distribution by 
virtue of inequality (1). 

Finally, observe that D takes values in [h*, oo). As in the proof of Theorem 
5 of [2], it follows from inequality (2) that V(X F ^) can be dominated by 
_Dp( 2 ), as required. □ 

The majority of Section 3.3 remains valid when the dominating process 
is constructed as above. The only issue is that by taking S = F in this new 
method we are no longer assured that S(h*) > m, where the set C* = {x : 
V(x) <h*} is m-small. Unfortunately, there no longer seems to be a simple 
way to ensure this since our attempts to increase S in the above always 
result in an increased value of h* . 

If it so happens that F(h*) > m for a given chain, then the original perfect 
simulation algorithm remains unchanged. If this is not the case, then the 
algorithm must be altered. It now becomes necessary, when Dq = h* , for D 
to dominate V(X) not at time o\ = F(h*) but at time 

a* = inf jo",- : On > ml. 

j>2 J J 
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This is an example of the composite nondeterministic sampling schemes we 
had originally hoped to avoid (cf. the comment before [1], Theorem 15]). 
Furthermore, we need to be able to couple target chains and dominating 
process at a* in such a way that the target chains may regenerate at this 
time (using the fact that C* is cr*-small). This unfortunately reduces the 
impact of the result, which is an issue that we are currently trying to resolve. 
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